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Abstract 

Nucleation, Growth, and Coarsening — A Global View on Aggregation 

by 

Joseph Farjoun 
Doctor of Philosophy in Mathematics 

University of California at Berkeley 
Professor John C. Neu, Chair 

We present a new model of homogeneous aggregation that contains the essential physical 
ideas of the classical predecessors, the Becker-Doring and Lifshitz-Slyovoz models. These 
classical models, which give different predictions, are asymptotic limits of the new model 
at small (BD) and large (LS) cluster sizes. Since the new theory is valid for large and 
small clusters, it allows for a complete description of the nucleation process; one that can 
predict the creation of super-critical clusters at the Zeldovich nucleation rate, and the 
diffusion limited growth of large clusters during coarsening. By retaining the physically 
valid ingredients from both models, we explain the seeming incompatibilities and arbitrary 
choices of the classical models. 

We solve the equations of our new model asymptotically in the small super- 
saturation limit. The solution exhibits three successive 'eras': nucleation, growth, and 
coarsening, each with its specific scales of time and cluster size. During the nucleation era, 
the bulk of the clusters arc formed by favorable fluctuations over a free energy barrier, ac- 
cording to the analysis by Zeldovich. The free energy barrier increases as more clusters are 
formed, and this signals the beginning of the growth era: no new clusters are created, and 
the expansion of the existing ones continues. The growth of the clusters slows down when 
the reservoir of monomers that fuels it is sufficiently depleted. This signals the onset of the 
final coarsening era. This is a competitive attrition process, of smaller clusters dissolving 
and fueling the further growth of the larger survivors. By resolving the preceding creation 
and growth eras, our analysis gives explicitly the characteristic time and cluster size of the 
coarsening era, and a unique selection of the long time, self-similar cluster size distribution. 
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Chapter 1 



A Global Model 



1.1 Introduction 

Nucleation refers to the aggregation of identical particles (monomers) into clus- 
ters. Its universality throughout physics, chemistry and biology is well known. References 
[KGT83j . |Kel91j . [NBOO] , [GWS+Ol] . [Is?9T] . |^pB02| , [LSHT] . |XH91j . |MG96j . |G^ON03j 



provide a lineup of the 'usual' (and some unusual) suspects. Also well known are the long- 
standing challenges that aggregation poses to modeling. Two classical models of aggregation 
due to Becker-Doring (BD) [BD35j . and Lipshitz-Slyozov (LS) |LS61j are incomplete and 
mutually inconsistent. 

In BD, clusters exchange particles with the surrounding monomer bath by a 'sur- 
face reaction', and it is assumed that the monomer bath around the clusters has uniform 
concentration. This is only possible with infinite diffusivity of monomers. While this de- 
scription is asymptotically accurate for sufficiently small clusters, the uptake of monomers 
by large clusters is strongly controlled by the diffusivity. LS describes cluster growth and 
shrinkage controlled by diffusion of monomers. In LS, the monomer concentration at the 
surface of a cluster is a prescribed function of the local curvature, generally different from 
the 'background' concentration, far from clusters. Hence, monomer concentration about a 
large cluster is nonuniform, and there is diffusive transport of monomer into or away from 
the cluster. This physics of LS leads to a prediction for cluster growth that disagrees with 
BD. Furthermore, LS is 'incomplete', in that it does not describe the initial creation of 
clusters from pure monomer. While it is generally accepted that BD is a model for small 
clusters, and LS for large, several questions remain. How to interpolate between the two 
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models? What is the characteristic size that separates 'large' from 'small'? What physics 
governs the growth in the intermediate scale? What globally valid model encompasses 
the whole evolution of clusters, from an initial state of pure monomer to the asymptotic 
self-similar distribution of large cluster sizes? 

The current chapter presents a new model that retains the essential physical in- 
gredients: the clusters gain and lose monomers by a surface reaction that depends on the 
cluster size and the monomer concentration seen on the surface. Monomers outside the clus- 
ter undergo diffusion with finite diffusivity. These ingredients give rise to a free boundary 
problem for the growth of a cluster that contains a new intrinsic cluster size, k^:, in addition 
to the well known critical size, kc- The critical size, kc, separates shrinking clusters {k < kc) 
from growing ones {k > kc)- The new cluster size, /c*, indicates the importance of diffusion: 
the new prediction for cluster growth asymptotes to BD for small clusters with A; <C A;* , and 
to LS for large clusters with k ^ k^,. In the former case, the diffusion effectively equates 
the surface density of monomer with the far-field density, thus, the surface reaction dictates 
the growth. In the latter case, growth is strongly limited by finite diffusivity. Furthermore, 
the new model of cluster growth interpolates between BD and LS for intermediate cluster 
sizes on the order of A;*. 

The smooth interpolation between BD and LS is crucial for a global model of 
aggregation that describes the whole process, from the initial creation of clusters from pure 
monomer, to the late stage growth-attrition process called coarsening. The essential idea 
is simple: if kc <C k^,, as expected in most cases, standard BD describes the nucleation of 
super-critical {k > kc) clusters and their growth while kc < k <^ k^,. The super-critical 
clusters rapidly grow to sizes k ^ k^:, and their subsequent careers are described by LS. 

Mathematically, we model this physics by a continuum approximation of the dis- 
crete kinetics. The continuum equations constitute a PDE signaling problem for the distri- 
bution r{k,t) of large {k <^ k^) clusters in the space of (continuous) cluster size k. At the 
lowest order of approximation, the cluster-size distribution satisfies an advection PDE, in 
which the growth rate {k vs. k) furnishes the advection velocity. The classical Zeldovich 
formula [Zel43j . which follows from BD, computes the creation rate of super-critical [k > kc) 
clusters. Since we assume kc "^k^, the Zeldovich formula gives rise to an effective source 
boundary condition on A; = 0. The initial state of pure monomers is expressed by a zero 
initial condition, r(A;, 0) = 0. Information about the amount of small (A; < kc) clusters is 
not expressed directly in r{k,t). Instead, using conservation of particles, we express the 
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amount of sub-critical clusters using an integral of r(k,t). 

Our theory does not handle nucleation that happens with kc on the order of k^,. For 
this we suspect that a new theory is needed, one that considers the discrete and fluctuating 
nature of the monomer bath, and does not resort to the diffusion equation, which arises 
from mean-field averaging. 

The chapter is structured as follows. In section [L2] we present a short summary of 
the classical microscopic aggregation theory (BD). We derive rate constants for attachment 
and dissociation of monomers from a cluster by using free energy and detailed balance 
arguments. The only difference from the classical theory is that it is based on the surface 
density of monomers, the density of monomers just outside the cluster, and does not assume 
that the monomer density is homogeneous. 

In section 3 we take into account the finite spatial diffusion of monomers. While 
still focusing on a single cluster, we connect the surface monomer density with a far-field 
monomer density, the nearly uniform concentration of monomers far from any cluster. This 
prescribes the growth rate of a cluster as a function of its size and the far-field monomer 
density. The standard assumption in diffusion limited aggregation is that the surface density 
corresponds to a critical cluster i.e., growing and shrinking are equally likely. This seems 
paradoxical for two reasons. First, the free energy of a cluster as a function of cluster size, 
k, has its global maximum at A; = kc- On the face of it, this seems to be an unstable 
equilibrium, but yet it is claimed that the cluster remains at the top of this equilibrium. In 
addition, if that is the value of the monomer density, how does the cluster grow or shrink? 

This paradox is another artifact of assuming uniform monomer concentration as 
in BD. It is deconstructed at the end of section 3, by an asymptotic analysis which exposes 
the stabilizing role of finite monomer diffusion: If monomer concentration at the surface 
of a large {k ^ A;*) cluster has large deviation from the critical value described above, the 
surface reactions rapidly absorb or expel monomers. Consequently, due to finite diffusivity, 
the surface monomer concentration undergoes a collateral adjustment towards the critical 
value and the rapid reactions are turned off. The growth rate due to diffusion is much 
slower and dictates the evolution of large clusters. 

In section [L^ we turn to the ensemble of all clusters. If the density of monomers is 
below a certain saturation value, an equilibrium exists, in which large clusters are extremely 
unlikely. For a 'super-saturated' ensemble, with monomer density exceeding the saturation 
value, clusters greater that a 'critical cluster size', kc, have a strong tendency to persist and 
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grow. In this super-saturated case, there is no equihbrium; the distribution of clusters is 
continuously changing. Initially, there is nucleation, which is the creation of super-critical 
clusters. The calculation of the nucleation rate based on BD is reviewed here. 

As stated before, this chapter proposes a PDE signaling problem for the distri- 
bution of cluster sizes that quantifies the complete evolution of the aggregation process. 
Section 5 contains the assembly of the signaling problem from the component parts in sec- 
tions 2-4. It has a peculiar nonlinearity, in which the advection velocity in the PDE and 
the boundary condition (BC) at A; = depend on the monomer density as a parameter. 
The monomer density can be written as an integral of the solution, and herein lies the 
nonlinearity. 

The nonlinearity makes the task of solving the equations difficult enough to warrant 
placing it in the next chapter. 

1.2 Classical Becker-Doring Model 

Becker-Doring theory (BD) imposes simplifying assumptions at the outset: The 
clusters are assumed to be uniformly distributed in a dilute 'bath' of monomers. This 
assumption is adjusted regarding the distribution of monomers in the next section, when 
we add diffusion. The clusters are assumed to change size only by losing or gaining one 
monomer at a time. Two large clusters do not fuse together nor does one cluster break into 
two. This can be justified hcuristically by noticing that the density of the large clusters 
is much smaller than the (already small) density of monomers, thus the probability of two 
large clusters interacting is small. In addition, the mobility of the large clusters is much 
smaller than that of the monomers, so they arc even less likely to stumble upon one another. 
Similarly, since the large clusters have a low mobility (relative to monomers) a cluster that 
breaks into two will, most likely, reconnect quickly, as the two parts remain close together. 

Another important assumption is that the only governing parameter of a cluster is 
its size. The shape of the cluster is assumed to be fixed. This assumption can be weakened 
to require that clusters of same size have the same binding energy and the same surface 
area. 

To derive the kinetic model of nucleation we introduce the essential quantitative 
ingredients: energy, free energy, and the rate constants of transitions between configurations. 
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1.2.1 Energy and Free Energy 
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Figure 1.1: A schematic cubic crystal sur- 
rounded by the monomer bath. 

The energy of cubic cluster with k monomers is 

ek^3e{k-k^. (1.1) 

Here £ is the binding energy of a single bond between two adjacent particles, the energy 
needed to break it. 

While the clusters are not assumed to be simple cubes, the general structure of 
the binding energy is expected to remain. Thus, we assume that there is a bulk energy 
constant, a > 0, and a surface energy constant, o" > 0, such that the energy of a A;-cluster 
is 

3 2 

Ek ~ kBT{ak - ^(rk3), for /c> 1, (1.2) 

where, ksT is the Boltzmann factor. The factor of | is added in hind-foresight, as it makes 
some formulas cleaner. Equation (jl.2p is only true asymptotically for large clusters. For 
small clusters, the separation between 'bulk' and 'surface' is artificial, and we do not expect 
(|1.2|) to be quantitatively accurate. In particular, ei = since the binding energy is the 
change in the energy from the unbound state and a cluster with one particle is unbound. 

Next we consider the free energy costs to create a fe-cluster from the monomer 
bath. The bath is characterized by the density pi of monomer, measured in units of i, 
where v is the volume of each monomer. In other words, pi is the volume fraction occupied 
by monomers 
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The free energy cost to create a fc-cluster from the monomer bath is 



9k = -^k - ksTk log pi 

« -ksT ({a + logpi)k - ^ak^ ) , for A: > 1. 



:i.3) 



Here, ksTlog pi is the chemical potential of a monomer in the bath. Rewriting the k^ 1 
asymptotic form of the free energy gives insight into the existence of a critical monomer 
density, ps = e~". We call ps the saturation density of monomers. Setting a = log ^ in 
(|1.3p gives 



, ^ / 3 2 pi 
gk ~ kBT -ak3 - fclog — 
V2 Ps 

Thus, when pi < ps the free energy increases with A;, allowing for an equilibrium. When 
pi > Ps, the free energy attains its maximum at the critical cluster size kc-, 



a 



, as Pi ^ 



We investigate the implications of this critical value later. 



1.2.2 Kinetics and Detailed Balance 




(1.4) 



Figure 1.2: The exchange of particles between a 
cluster and the monomer bath. piCk is the rate 
at which a monomer gets added to the cluster, 
and dk is the rate at which monomers leave. 



To study the rate of change in the cluster size we need to model the allowed 
reactions. BD allows two types of reactions:(l) A monomer in the bath can join the cluster; 
(2) A particle on the surface of the cluster can dissociate from it and enter the bath. 
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Mathematically, the two reactions are modeled as independent Poisson processes. Let us 
call the rate at which particles leave a (A; + l)-cluster d^, and the rate at which particles 
join a A;-cluster CkPi- We work under the assumption that the solution is dilute, that is, 
pi <^ 1. It is reasonable to expect that in the dilute limit, the growth rate is proportional 
to pi. Thus, we include pi in the growth rate, CkPi, so that both Ck and dk are independent 
of pi. 

Initially we assume a uniform pi, so it is a global parameter. When we add 
diffusion, we are a little more careful and take pi to be the density at the surface of the 
cluster. 

Since the dissociation happens on the boundary of the cluster, we expect to be 
proportional to the surface area, which, in turn, is proportional to k"^^^, 

2 

dk = ujk^ . (1.5) 

Here, uj includes a per-particle dissociation rate and a geometric factor. 

A standard detailed balance argument relates Ck to dk- Let pi be the value of 
monomer density so that a /c-cluster is in equilibrium with the monomer bath. That is, the 
cluster has no net tendency to grow nor shrink, hence, the adsorption rate, piCk, should 
match the emission rate, dk- The physical requirement for the equilibrium is that the free 
energy is unchanged when a monomer is taken from the bath and added to the cluster, 

gk+i - 9k = -£k+i +£k- ksT log pi = 0, 

hence, 

Pi = e ''B^ . 

The detailed balance relation between Cfc and dk is therefore, 

Ck = e ''bt dk- (1.6) 

We use the adsorption and emission rates, CkPi and dk, to derive a kinetic equation for the 
expected change in size of a cluster. In a small time span 5t, the expected change in cluster 
size, 6k, is: 

{6k) = {ckPi - dk)6t- (1.7) 

Using the model for the binding energy of large clusters (|1.2p and relation (|1.6p between Ck 
and dk, equation (jl.7p can be rewritten as 

{6k) w UJ ( ^^~'^' kl - ak^i^ 6t, for » 1. 
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In the appendix we show that if the variance in k is much smaller than {k) initially, it will 
remain so, as long as {k) is bounded away from the critical size kc in (jl.4p . In this case we 
approximate the evolution of k{t) for a given cluster as deterministic, and governed by the 



Here, r/ is the super- saturation, defined by 



The super-saturation in ()1.9|) is generally a function of time, due to the exchange of particles 
between clusters and the monomer bath. By conservation of the total particle density, the 
average value of pi, and hence 77, is determined from the densities of all clusters with k > 2. 
Within the framework of BD, which assumes that the monomer density is uniform, we 
simply set rj to this average value, and, in this sense, (jl.8p is the BD prediction for cluster 
growth. However, if the diffusivity of monomers is finite, the density of monomers seen at 
the surface of the cluster will be different from the average value far away. We propose that 
p.Sp holds generally, with r] equal to the value of super-saturation seen at the surface of the 
cluster. Equations (jl.Sp and (|1.9p expose the kinetic significance of the saturation density 
Ps = e~", and the critical cluster size kc in (jl.4p . which we rewrite using r], 



If the surface value of pi is less than ps, all the clusters shrink, regardless of k. For pi > ps, 
i.e. 7] > 0, the critical size, kc, separates growing, super- critical clusters {k > kc) from 
shrinking, sub-critical ones {k < kc). 

It should be noted that while the expected change in size of sub-critical clusters 
is negative, there is a small probability for a sub-critical cluster to grow and become super- 
critical. The rate at which this happens is estimated by the Zeldovich formula, to which we 
come back in section 11.41 First we add finite diffusion of monomers to the model and see 
how it affects the growth rate. 



ODE 




(1.8) 




(1.10) 



1.3 Adding Monomer Diffusion 



Diffusion is the usual model of transport for monomers in the bath. For a finite 
diffusion coefficient, D, the exchange of particles between a cluster and the bath directly 
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Figure 1.3: The cluster is surrounded by a 
inhomogeneous monomer bath. The flux of 
monomers in the bath is diffusive, and the clus- 
ter reacts to the local monomer density. 



outside of it leads to non-uniform monomer density. Hence, the monomer concentration, 
pi, is a function of position and time that satisfies the diffusion PDE. 

To add diffusion to the BD model, we make a couple of additional assumptions. 
The cluster is taken to be spherical, filled with monomers, each taking volume i'. Thus, the 
number of particles in the cluster, k, and the radius of the cluster, a, are related by 

Air 

ku=—a^ (1.11) 
3 ^ ^ 

The monomer density is assumed to be radially symmetric (with the cluster centered at the 
origin.) Thus, the density, pi{r,t), satisfies the radially symmetric diffusion PDE in IR'^ 

dtpi = ^drir^drpi) inr > a, (1.12) 

There are two BC at r = a: one is the kinetics due to BD as in equation (jl.Sp . 
That is, pi{a,t) is related to k and k by 



k = uj \ rj{a, t)k'A — fjfca J , (1-13) 

where, 

pi{a,t)-ps 
7]{a,t) = - 



Ps 

is the super-saturation seen at the surface of the cluster. 

The second BC results from conservation of particles. The particles that are added 
to the cluster can come from two possible sources: particles in the solution surrounding the 
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cluster, which join the cluster as it engulfs them, and particles added by the diffusive flux. 
This can be put in a simple equation: 

ail - pi{a,t)) = D{drpi){a,t) = Dps{drr])ia,t). (1.14) 

Finally, there is a BC at oo: the monomer concentration has the asymptotically 
uniform value Poo far from the clusters. 

Pi — > Poo as r ^ oo. (1-15) 

In the context of the full aggregation problem, poo is generally a function of time, which fol- 
lows from overall conservation of particles. Equations ()l.llfffTT5]) constitute a free boundary 
problem (FBP) for a{t) and pi{r,t) in r > a. 

The analysis of the FBP begins by identifying suitable non-dimensional variables. 
We assume that the local super-saturation 

r]{r,t) = . 

Ps 

is uniformly small in r > a. Hence, we introduce e as a gauge parameter for r]{r, t) and we 
replace 7/(r, t) in (jl.llffLTS]) by er/(r, t). The analysis of the FBP is based on an e ^ limit 
process. However, we do not yet take the e ^ limit. It remains to determine the scaling 
of the other variables k, r, a and t with e. These follow from simple physical balances: in the 
'kinetic' BC (|1.13p . r]{r,t) has order of magnitude e, so the two terms of the RHS balance 
when k has magnitude (f • Notice that this is the critical cluster size (jl.lOp for jy = e. 
The scalings of ry and k are recorded in the scaling table 

Scaling Table 
Variable r] k a,r t 



Unit e (f)^ 



The unit of r and a is the radius of a cluster with k = (f )^ particles. The unit of t is chosen 
to balance the LHS of the BC (jl.lSp with the two terms on the RHS. 
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The equations of the non-dimensional FBP that follow from ()1.11HLT5]) are 

47r o , , 

k = —a^ (1.16) 

dtV = — ^ — \drir^drr]), in r > a, (1.17) 
k = r]{a,t)k^/^ - k^/^ , (1.18) 

T] — > r/oo as r ^ oo. (1-20) 



d (l - p.(l + er/(a, t))) = -^^-^{drr]){a, t) (1.19) 



Here, r/oo is the asymptotically uniform value of ri{r,t) as r ^ oo. The characteristic time 
of its variation is assumed to be comparable to or larger than the unit of time from the 
scaling table. 

The dimensionless constant ^^/s can be interpreted as a ratio of characteristic 
times for two different physical processes. Recall that ^ is the characteristic time for a 

2/3 

monomer on the surface of a cluster to dissociate into the solution. The ratio ^^-p- is the 
characteristic time for a monomer to diffuse a distance comparable to its own size. The 
conventional assumption is that the "dissociation time" is much longer than the "diffusion 
time," so that 

1 ^ 

In this limit, the diffusion equation (jl.l7p reduces to a radial Laplace equation 

drir'^drf]) = 0. 

The solutions with tj = r]oo as r — > oo are 

r]{r,t) = r]oo + (via, t) - Voo)-- (1-21) 

r 

The time dependence of r] is implicit due to the time-dependence of its values at r = a and 
r = oo. 

In the dilute limit, with ^ 1, the conservation equation (|1.19p reduces to 

a = -^^^ 
Substituting ry from (|1.2ip . this is becomes 



■{drri){a,t). 
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Equation ()1.22p can be converted into an equation for k using (jl.l6p . This gives 

k = eiJ.{r]oo -r]{a,t))k^, (1.23) 
where is the dimensionless constant: 



fi={3- levr^)^ 



D \ (P. 



K = {anf = (3 • Wtt^)^^. (1.26) 



Since /i is a product of a large number, ^^2/3 ; ^^^d a small number, it can be large or 
small. We therefore entertain any value of fi. 

The two equations for k, (jl.lSp and (jl.23p . involve the super-saturation at the 
surface of the cluster, ri{a,t). Solving for k and rj{a,t) gives 

^(a,t) = ^^. (1.25) 

ODE p.24p indicates a second characteristic cluster size besides kc- In units of (f )^ this 
cluster size is 

k^ = {efif. 

In the original variables, k^, is a combination of basic physical constants: 

\3 _ i«_2-. 

Notice that for k <^ k^ equation (jl.24p asymptotes to BD. For k ^ k^, the asymptotic form 
of (fOil) is 

k-efi (^r]ook^i - l) . (1.27) 

Restoring original units, (I1.27P becomes 

k = d (r]ookri - a) , d = (3 • lOvr^)! (1.28) 

which is the standard result for diffusion limited growth (DLG) [LS61] . Equation (jl.25p 
shows how the surface value of super-saturation differs from the uniform value, r]oo, far from 
the cluster. Notice that it is a function of k. We convert it into an equation for Sk, the 
value of monomer density seen at the surface of a A;-cluster (again, in original units): 

This will be important, when we examine the whole ensemble of clusters, and formulate 
evolution equations for cluster densities. 
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1.3.1 Physical Meaning of 

We show that k^, is the characteristic size of clusters for which finite diffusion 
induces a significant relative difference between rjoo and rj{a, t). That is, 67] = rjoo —Tj{a, t) is 
comparable in magnitude to r^oo- A simple examination of two physical balances is sufficient. 
First, the cluster's growth rate balances the diffusive infiux of monomers. This is expressed 
by 

k = a^D^^l^. (1.30) 

a 

Here, the equality means 'order of magnitude balance'. In the RHS, ^6r] is the difference 
between monomer densities at 00, and on the surface, expressed in the conventional unit 
of 1/volume. For quasi-static diffusion, the diffusion zone about the cluster of radius a has 
thickness a, so -^-^ — estimates {drPi){a,t) and the influx of monomers per unit area into 
cluster is estimated by multiplying this by D. Finally, multiplying by the area, proportional 
to a^, gives the cluster growth rate, k. 

Second, the magnitude of k as dictated by the surface reactions (jl.l3p is 

k = ojr]kl. (1.31) 

Enforcing the equivalence of (11.30p and (11.3ip and using k = va^ (order of magnitude 
equality again), we find 

5-n ujv^l'^ , 1 
— = — — A;3. 

And we see that 5ry is comparable to r\ when k is comparable to /c* = {crp)^. 

1.3.2 Critique of DLG and its 'paradox' 

We briefiy examine the 'traditional' derivation of ODE (jl.27p for DLG, within the 
framework of the non-dimensional free boundary problem (jl.l6HL20|) . Given rj{a,t), equa- 
tion (jl.23p gives the growth rate of the cluster that follows from diffusive flux of monomers. 
In the traditional analysis of DLG, r/(a, t) is chosen so that the cluster is in equilibrium 
with the monomer bath that surrounds it. Under the current non-dimensionalization, this 
'critical nucleus' BC reads 

rj{a,t) = k~^. (1.32) 
Substituting (fL32]) for r/(a,t) in (fr23]l leads directly to the ODE (fOTl) . 
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By inspection, we see that the traditional equations (jl.27l ll.32p arise by taking 

the e/i ^ limit of (jl.24l I1.25P with k fixed. The alternative limit process, ^ > oo with 

= Efj, fixed is more physical: the value of efi is set by material properties and initial 
conditions, and we expect that clusters eventually grow to sizes k ^ k^:. We have already 
seen that the ODE (jl.24p for k converges to the DLG result in this limit, but the expression 
(jl.25p for the surface value of monomer density does not converge to the DLG boundary 
condition ()1.32p . Instead, 

r]{a,t)^{l + Hr]^)k~^, (1.33) 

which has an additive term ^rjoo in the prefactor of k^^^^ not present in (jl.32p . A mathe- 
matical critique of the 'critical nucleus' boundary condition (|1.32p is simple; it results from 
formally neglecting k in the LHS of the 'surface kinetics' boundary condition (jl.lSp . In our 
result, k balances the RHS, even in the limit — > oo, resulting in (11.33p . 

Recall that the 'critical nucleus' boundary condition in traditional DLG looks para- 
doxical because the 'cluster sits on top of a free energy maximum'. A 'lazy' deconstruction 
might say: "Nothing to explain, the critical nucleus boundary condition is simply incorrect 
in the (more physical) limit ^ oo with k^ fixed." Another easy explanation looks at 
the free energy. The free energy (jl.3p refers to a simple cluster surrounded monomers of 
uniform density, whereas the actual kinetics we consider involves a non-uniform density 
pi{r,t) in r > a. The actual free energy takes into account the functional dependence of 
pi{r,t) in r > o. These remarks indicate that the 'paradox' in its original form is nai've. 
Nevertheless, it points to some physics that is not expressed in the quasi-static model (|1.24l 
I1.25P of cluster growth as it stands. 

Suppose that we place a cluster of size k into a uniform monomer bath that has 
the 'wrong' monomer density, not equal to the surface value Sk given in (|1.29p . In order 
for our model to be plausible, the surface value of pi should rapidly relax to Sk in (|1.29p . 
We now show that the full free boundary problem (|1.16I[L20|) implies such a relaxation 
transient. 

1.3.3 The Stability of the Free Boundary Problem 

The relaxation transient is characterized by a balance of the time and space deriva- 
tives in the diffusion PDE (jl.l2p . Hence the characteristic time of the relaxation transient 
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is 

tr ^ , (1.34) 

where a is the cluster radius. The relative change of the cluster radius in this characteristic 
time is small: from (11.140 . d has the order of magnitude ^^"^ . Hence, the relative change 
of radius in time tr has magnitude eps- The small relative change in cluster radius means 
that the cluster radius is asymptotically constant during the relaxation transient, and it 
remains to derive from the full free boundary problem (jl.llffLTS]) a reduced boundary 
value problem for r/(r, t) in r > a, with a fixed. We use the previous units in the scaling 
table for all variables except time t. For t we use tr in (11.340 with a replaced by the 
characteristic cluster radius, ji^^^^- The reduced boundary value problem is 

dtT] = —dr{r'^dr'q) in r > a, (1.35) 
\{drr]){a,t) = r]{a,t) - k-'^i , (1.36) 
ry(r, t) ^ r?oo as r ^ oo, (1-37) 

in the limit e ^ 0, and A = (^j ^ fixed. The far- field super-saturation, ryoo; is assumed 
to vary on a characteristic time much longer than t^j so "q^o is effectively constant. 

The time-independent solution of (jl.35ffL37|) for 7/(r, t) is (jl.2ip with rj on r = a 
given by (|1.25p . We show that this time-independent solution is asymptotically stable. We 
notice that 

A a? / 1 \ 2 

^^2] ^^i9rvfdr + -{rj{a,t)-k-^j (1.38) 

is a Lyapunov functional for equations ()1.35ffL37|) . The time derivative of E, 

l-CO 

E = —A / r^?7j dr, 

J a 

is found by time-differentiation of (|1.38p . integration by parts, and use of the PDE (jl.35p 
and BC (11.360 . Since E is positive definite, and ii^ < 0, it follows that ri{r,t) converges to 
the time-independent solution. We conclude that if the surface monomer concentration is 
initially different from the quasi-static value ()1.25p . it relaxes to it in characteristic time tr- 

1.4 Evolving Distribution of Cluster Sizes 

The kinetics equation (11.240 requires tjoo, the super-saturation far from any cluster. 
To find rjoo we look at the joint evolution of all the clusters, each assumed to follow the 
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dynamics in (jl.24p . The clusters are coupled by the combined effect they have on the 
monomer density, and consequently, on the super-saturation. 

Let pkit) be the average spatial density (in units of ^) of A;— clusters at time t. 
Assuming that the total particle density has a fixed value, p, the pk satisfy a particle 
conservation equation. Since a fc-cluster is made of k particles, the total particle density, p, 
must satisfy 

oo 

p = Y,kpk- (1.39) 

k=l 

With the help of (jl.39p . the space averaged super-saturation can be rewritten as a function 
of the cluster densities pk with k >2: 



Pi- Ps P- Ps 



^ oo 

-Y^kpk. (1.40) 



Ps Ps p. 

In the dilute limit with inter-cluster distances much greater than cluster radii, we expect 
that the super-saturation is asymptotically uniform, throughout most of the monomer bath 
far from clusters. In this case, that asymptotically uniform value, r]oo{t), should be well 
approximated by the spatial average (jl.40p . 



oo 

P- Ps i 



^ oo 

-Y^kpk. (1.41) 



We turn to the evolution of the densities. The pk obey kinetic equations associated 
with the reactions 

/c-cluster + monomer <^ (A; + l)-cluster. 

The equations are 

Pk=jk-i-jk, (1-42) 

for k > 2, where the discrete flux jk is the net rate of creation of a {k + l)-cluster from a 
fc-cluster, 

jk = CkSkPk - dkPk+i- (1-43) 

As before, dk is the rate constant for shedding a monomer from the surface of a (A; -|- 1)- 
cluster. For A; ^ 1 it has the asymptotic behavior (|1.5p . proportional to surface area. The 
prefactor CkSk of pk in (jl.43p is the rate constant for adding a monomer. Recall that Sk 
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is the value of monomer density seen at the surface of a /c-cluster, given by (jl.29p . and 
Cfc is related to by the detailed balance condition (jl.6p . A more explicit formula for 
displaying the /c-dependence of surface monomer concentration and detailed balance is 

jk = dk{e SkPk - Pk+i ■ (1-44) 



Equations ()1.421 ll.44p can be summarized as discrete advection-diffusion equations, 



Pk + D- 



4(1-6 ''bt sk] Pk- dkD^pk 



for /c > 2. (1.45) 



Here D^,D^ are, respectively, the forward and backward difference operators. In these 
equations, the surface monomer density Sk contains the super-saturation ryoo as a parameter, 
and rjoo is connected to the pk according to (|1.4ip . So we see explicitly how the cluster 
densities are coupled to each other via the super-saturation. The /c-dependence of Sk induced 
by the finite diffusivity of monomers is the essential difference from classical BD. We recover 
classical BD by taking p ^ oo, which in turn results from ^^/^ — > oo. Then Sk in ()1.29p 
reduces to ps{l + rjoo), which is the uniform value of pi assumed in classical BD. 

1.4.1 Equilibrium 

Equilibria are time independent densities, pk, so that all the fluxes jk are zero, 
and the sum (jl.39p . which gives the total particle density, is convergent. Setting jk = in 
(jl.44p gives a recursion relation that determines pk from pi, 

p^. = p^e^, for k>2. (1.46) 

Here, we used Sk = pi for all k, since pi should be uniform in the equilibrium case. Substi- 
tuting these Pk into ()1.39p gives, 

P = p{pi)^Y.kp',e''BT. (1.47) 

k=l 

Hence, equilibria exist for monomer densities pi so that this series converges. For large 
values of k, the binding energy can be written as 



Ek ~ keT (^ak - ^ctA:^/^^ 



Therefore, convergence happens for 
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In other words, equilibria exist only if the super-saturation is non-positive. The largest 
value of total particle density for which there is equilibrium is obtained by setting pi = ps 
in (|1.47p . We denote this critical particle density by pc, 

Pc = p{Ps). (1.48) 
Since ps <^1, the first few terms of the series give a close approximation of Pc- 

1.4.2 Zeldovich Nucleation Rate 

For positive super-saturation r], there is the critical cluster size k = k^, for which 
the free energy cost to assemble a /c-cluster from dissociated monomers is maximized. In 
the small super-saturation limit r/ ^ O"*", and k ^ 1, it follows from (|1.3p and (|1.9p . that 

gk ~ keT (^^akl - ai]^ , (1.49) 

and that the free energy cost (in units of /c^T) of the critical cluster is asymptotic to 

k zri'^ 

For small super-saturation, the free energy cost is high, and an analogy with the famous 
Arrhenius rate suggests that super-critical nuclei with k > kc are produced at a rate pro- 
portional to the exponential e~^. Since this proposed creation rate is exponentially small as 
rj ^ 0, one might expect that after some initial transient, quasi-static but non-equilibrium 
values of pk are established for k on the order of kc, in which the discrete fluxes jk in (jl.43p 
are asymptotically equal to a uniform value j. This j is the creation rate of super-critical 
nuclei, proportional to e~^. These essential ideas of nucleation kinetics are set forth in a 
famous paper of Zeldovich, on the nucleation of vapor bubbles for under pressurized liq- 
uid [Zel43] . His starting point is a discrete system of kinetic ODE's like BD, but he first 
passes to a PDE limit of the ODE's and calculates the nucleation rate from the PDF. Here, 
we implement the essential Zeldovich ideas, but within the framework of the discrete BD 
ODE's. 

We work in the limit k^ <C fc*, so for k on the order of kc there is negligible 
difference between the surface value, s^, of monomer density in (jl.29p . and the uniform 
value, pi, far from clusters. We show this: in p.29p for s^, we see that Sk ~ Psi^ + Voo) = Pi 
if J < f?oo and ^ < 1. For kc w ((t/??oo)^ in (fTTOl) and k^ = {apf in (fr26D . we find 
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~ [fi-q )^ • K ^ implies ^ <C r/oo- Furthermore, for /c on the order of kc, is on 
the order of ^ < 1. 

In (jl.44p . we replace Sk by pi, and by j, to obtain a recursion equation that 
determines the pk for k >2 from pi. We write it as 

pt_pm = ^M. (1.50) 

Here, is the free energy cost of a /c-cluster, given in (|1.49|) . and denotes the solution 
()1.46p of the homogeneous recursion relation with j = 0, and pi = pi. For positive super- 
saturation r/, — > oo as A: ^ CO, and we expect that the pk in (jl.SOp have ^ ^ as 
/c — > oo. Summing (jl.SOp over k gives a formula for j. On the LHS, we get a telescoping 
sum with value 

^_ ^ = 1-0 = 1. 

Pi fc-*oo Pk 

Hence, 



El gfc 
e'=sT. (1.51) 
dk-i 

k=2 ^ ^ 

In the RHS, decreases linearly with k as k ^ oo, so the series converges. In the small 
super-saturation limit — > 0^, we expect that the sum on the RHS is dominated by terms 
with k near k^ ^ {-] , where g^ attains its maximum. The relevant approximation to g^ 



as ?7 ^ 0"*" and k is on the order of kc is given by ()1.49p . Also, dk ~ ujk"^^'^ as in (jl.Sp . 
Hence, (jl.Sip has the asymptotic approximation 

• °° 2 

1 = ^^fc-iei"'^^-^^. (1.52) 

^ k=2 

The final step is the approximation of the sum by an integral, and evaluation of the r/ — > 0"*" 
limit by the saddle point method. This leads to the approximation of j, 

j lua — e (1.53) 
V Qtt 

1.5 Advection Signaling Problem 

We examine the aggregation process, starting from a super-critical density of par- 
ticles, p > pc all in the form of monomers at time t = 0. That is. 



pi(0) = p>pc, Pfc(O) = for > 2. 



20 




Pi(0) 



Pc 




> 



ps ps{l + r]*) Pi 



Figure 1.4: The darkened curve is the graph of p{pi) in 
< pi < Ps- The dashed hne is its hnear interpolation into 
Pa > rhog. It intersects p = pi(0) > pc at pi = ps(l + 
The arrow labeled i represents the decrease of pi during 
ignition, and the arrow a represents the decrease during the 
subsequent aggregation process. 



There is an initial transient, called ignition, in which the first supercritical clusters appear. 
A detailed analysis of the ignition transient appears in a paper by Bonilla et al. |NBC05) . 
Here we give a brief summary. First, sub-critical {k < kc) clusters are created with quasi- 
static densities close to the values pk in (11.46p . Of course, the value of pi becomes less 
than pi(0) = p, since these sub-critical clusters are created from monomers. Hence, the 
appearance of the sub-critical quasi-static densities is accomplished by the decrease of super- 
saturation from an initial value of ??(0) = to a smaller value, which we denote by rj^,. 
In appendix IA.21 we show that rj^, is related to p — pc by 



as p ^ p^. Here, p = p(pi) is the equilibrium relation p.47p between pi and p for 
< pi < Ps. Figure [131 is a visualization of the relation ()1.54p between 77=,, and p — Pc- The 
'completion' of the quasi-static densities for 1 < k < k^ = is accompanied by the 

appearance of the first super-critical clusters with k > kc- The rate of creation rises from 
zero to the Zeldovich rate in (jl.53p with rj = rj^,. We assume rj^, is so that kc <^ k^, so the 
Zeldovich rate indeed applies. 



Now our focus shifts to the evolving distribution of the super- critical clusters. The 



p- Pc^ 'n*PsP'{Ps) 



(1.54) 
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model we present here has three physical ingredients. Two of them are the growth of the 
clusters, and their creation. Both processes contain the super-saturation as a parameter, 
and the remaining ingredient is the connection of super-saturation to the distribution of 
cluster sizes. 

First, growth. We expect a predominance of super-critical clusters with /c ^ fc* 
that undergo diffusion limited growth. Here is the heuristic argument for not resolving 
size scales comparable to k^, or smaller: once a cluster achieves super-critical size, with 
^ (^^^ appendix lA.ip it continues to grow nearly deterministically, as shown 
in section 3. Since the Zeldovich rate is exponentially small in r], an exponentially long 
time elapses before the super-saturation shows an significant decrease below the effective 
initial value, 77*, established during ignition. In this exponentially long time, we expect the 
super-critical clusters to grow to sizes k ^ k^,, the regime of diffusion limited growth. 

We assume that the characteristic cluster size k corresponding to significant vari- 
ations of the densities pk for /c ^ /c* ^ 1 is itself much larger than A;* , and this motivates 
a continuum limit, 

Pkit)^rik,t). (1.55) 

Here, r{k, t) is a smooth function of its arguments, with the characteristic scale k of k much 
larger than A;*. Substituting p.55|) for into the discrete advection-diffusion equations 
()1.45p and using the assumed largeness of k, it follows that r{k,t) asymptotically satisfies 
the advection PDE, 

dtr + dk{ur)=0. (1.56) 

Here, the advection velocity u = u{k, rj) is identified from ODE (jl.28p for diffusion limited 
growth. We have, 

u{k,r]) = d(^r]k^ - , (1.57) 

where t] = ri{t) is the 'background' super-saturation, far from any cluster. 

Next, creation. In the analysis according to Zeldovich, recall that the discrete 
flux jk in (jl.44p is asymptotically uniform for k on the order of kc, with value j given by 
(jl.53p . Here, we make the stronger assumption that the range of k with the asymptotically 
uniform value j of jk extends to a scale of k much larger than A;* but smaller than the 
characteristic cluster size k associated with the continuum limit (jl.55p . In this case, we 
expect an asymptotic matching between the continuum limit of jfc, given by 

jk-u{k,ri)r{k,t), (1.58) 
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and the uniform value, j, of jk, in some overlap domain of cluster sizes k much larger than 
A;=K, but much smaller than k. Since k ^ k^, ^ k^ = (^-^^ in the overlap domain, the 
dominant component of u in (I1.57P is drjk^/^ and (jl.SSp reduces to 

jk ~ drjk^/^r{k,t). 

in the overlap domain. Hence, we propose the effective boundary condition on r{k,t) on 
k = are 

I 3 

drjk^/^ r{k,t) ^ i = ujj —e"^ as k ^ 0. (1.59) 

V 67r 

If 7/ = ri{t) is known, the advection PDE (11.56P and creation BC (jl.59j) lead to a simple 
determination of r{k,t), starting from the initial condition of pure monomer, r{k,0) = in 
A; > 0. 

It remains to connect ri{t) to r{k,t). In the particle conservation identity (|1.39p . 
we now distinguish between sub-critical and super-critical clusters sizes. That is, 

l<k<kc kc<k 

The sub-critical sum is approximated by substituting equilibrium values ()1.46[) for pk based 
on super-saturation r], 

Pk^Pk = {l + v)^Pse'^, 

and then taking the limit ry — > 0. The details here are a re-run of the calculation in appendix 
1X21 We get 

J2^Pk ^ Pc + VPsp{ps)- 
k=l 

The super-critical sum is approximated by the integral, 

r oo 

/ kr{k,t)dk. 
Jo 

Here, the lower limit is /c = and not kc, because the characteristic scale, k, of k in r(k,t) 
is much larger than kc- In summary, the conservation identity (jl.47p takes the asymptotic 
form 

/• CO 

p ^ Pc + VPsP'iPs) + / kr{k,t)dk. (1.60) 
JO 

Notice that if r = 0, which corresponds to negligible super-critical clusters, (jl.GOp reduces 
to ()1.54p with 7] = rji,. Hence, r/(t) has the effective initial condition r/(0) = 1]^,. 
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Thus, the signahng problem for r{k,t) consists of the advection PDE (jl.56p . the 
creation BC ()1.59p . and the functional dependence of the super-saturation, r/ on r{k,t) in 
(ll.60p . This signaling problem is nonlinear because rj, as a functional of r, appears in the 
advection PDE and, in addition, in the exponential creation rate in the BC. 

1.6 Conclusions 

The discrepancy between the two accepted model of nucleation — the surface reac- 
tion model derived by Becker and Doring, and the diffusion limited growth model due to 
Lifshitz and Slyozov — has now been resolved. Our new model, in which clusters interact 
with diffusing monomers by a surface reaction, predicts both models as limit cases. Al- 
though it provides a growth rate for clusters of all sizes, from small to large, its limiting 
behaviors are of special interest for us. In the limit of small clusters, the BD model emerges, 
and we can derive the Zeldovich creation rate of super-critical clusters using the BD kinet- 
ics. In the limit of large clusters, the diffusion limited growth model emerges, from which 
we derive the same evolution PDE for the density function, r(k,t), as LS. 

Not all our findings corroborate the classical assumptions and results. We find 
that the 'common wisdom' about the monomer density at the cluster surface is wrong. In 
the classical DLG model, the surface density is chosen so that the surface reaction is in 
equilibrium. In our model, the balance between the surface reaction and diffusion leads to a 
surface monomer density substantially different from the 'equilibrium' value. Despite this, 
the new model still predicts the same growth rate as the classical DLG model (for large 
clusters) . 

Another piece of 'common wisdom' is that BD pertains to small clusters and DLG 
to large ones. In our model, there is a new characteristic cluster size, k^,, that separates 
small (/c <C A;*) from large {k ^> k^,). This is useful in providing definite predictions for the 
validity of the model. For example, since nucleation happens around the critical size kc, 
and the Zeldovich creation rate assumes the small cluster limit of the model, the validity of 
the Zeldovich formula requires kc <C k^:. The case where kc is of the same order or larger 
than A;* is not covered by the current thesis nor by classical Zeldovich nucleation theory. A 
separate investigation is required for this case, one that goes beyond or treatment of the 
monomer bath as a smooth 'mean field' and really accounts for its discrete and fluctuation 
nature. 
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The monomer density around the cluster satisfies the diffusion equation and a 
mixed boundary condition at the surface of the cluster. To calculate the flux of the 
monomers into or out of the cluster, our model assumes that the monomer density in 
the vicinity of the cluster is a quasi-static solution of the diffusion equation. We show that 
the quasi-static 'surface' value of monomer concentration in our model is stable: starting 
from general initial conditions, the concentration of the monomer bath rapidly relaxes to 
the quasi-static approximation. 

The preceding insights are codified in a signaling problem for the cluster size distri- 
bution. It consists of an advection PDE consistent with diffusion limited growth, and a BC 
consistent with the Zeldovich creation rate. The PDE and BC contain the super-saturation 
as a parameter. The super-saturation, in turn, is a function of the cluster size distribution 
as dictated be the overall conservation of particles. Hence, the full signaling problem is 
non-linear. In particular, the Zeldovich creation rate is exponentially small as the super- 
saturation goes to zero. Therefore, accounting for small corrections to the monomer density 
is important, especially in the initial stage when new clusters are being created. 

Given an initial condition of pure monomers, quasi-static densities are formed 
during the ignition phase [NBCOSj . whereby small, sub-critical clusters are created from 
the monomers. The ignition transient is a precursor to the nucleation of super-critical 
clusters. Therefore, once the nucleation 'has ignited', the monomers density is lower than 
the its original value, prior to the ignition transient. The signaling problem starts from 
an effective IC that explicitly accounts for the loss of monomers during ignition. This is 
especially relevant in our case, where the monomer-density is close to the saturation density, 
Ps, and the dependence of the PDE and BC on the monomer density is very strong. 

Such is the state of the theory. We now conclude the conclusion by a small trespass 
into the domain of accountability: do the material parameters of real aggregation processes 
cooperate with the various assumptions of the modeling? We start with the characteristic 
cluster size k^,, most conspicuously present in our model. We examine physical constants 
associated with an aqueous solution of calcium carbonate, CaCOs. Its solubility is rather 
small, so the basic requirement of 'diluteness' is satisfied. Admittedly, CaCOs dissolves into 
positive (Ca^) and negative (CO3 ) ions, which is not reflected in our aggregation model 
with identical particles. But we are examining crude order of magnitude estimates, so will 
not be distracted by the inconvenience, and we formally consider a positive-negative ion 
pair as a 'monomer.' 
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Formula (jl.26p for fc* contains molecular volume the saturation density of 
monomer ps, the diffusion coefficient D of monomer in solution, and the dissociation rate 
constant oj. Estimates of v-,pc, and D are readily found in a chemical handbook |Lid05j . 
We use pc as an order of magnitude estimate of ps- The dissociation rate, w, is much more 
elusive. Here, we indulge in the activation energy model similar to Kelton's, in his review 
of glass-crystal transitions |Kel91] . The model is summarized by the formula 

o. = ^e-^ (1.61) 

2/3 

Here, ^^^is the characteristic time for a monomer to diffuse a distance comparable to its 
own size, and (3 is an 'activation energy of dissociation,' in units of ksT. Inserting (jl.6ip 
for uj into (jl.26p . we find 

3 



or, using ps = e 



fc, = (3- 167r2)e3(^-"). (1.62) 

Formula (jl.6ip has other applications for us, besides the estimate (|1.62p of k^:. For 
instance, recall that the quasi-static limit of the FBP ()1.2ip is based on the 'diffusion time' 
much smaller than the dissociation time ^, so we required ^^/3 ^ 1- From (|1.6ip we 

get 

^ (1.63) 



and so an activation energy kBT(3 on the order of a few KbT is sufficient. 

Next, we examine the criterion ^ ^ 1 for the validity of the Zeldovich nucleation 
rate. Using (jl.lOl [L26p . we find this implies a bound on the super-saturation. 

Using the crude 'cube' model of bonding energy (jl.ip we estimate a by, a |a so the 
criterion on the supersaturation becomes 

7/>ae"-^. (1.64) 



A lower bound on the super-saturation might seem like a problem, as we expect 
to asymptote to zero in the long-time limit of an aggregation process. However, significant 
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nucleation occurs at an early and relatively brief phase, so (jl.64p should apply for the initial 
super-saturation. In late stage coarsening, rj is much smaller than the RHS of (|1.64p . but 
significant nucleation is not happening then. 

Numerical evaluation of k^,, or the RHS of (jl.64p require actual values of a and 
p. It is relatively easy to find a, as we have a = — log ps and ps is essentially the volume 
fraction of monomer in saturated solution. For instance, from the data in |Lid05j . we find 
that 4.1 X lO^^cm'^ of solid CaCOa is soluble in lOOcm^ of water at room temperature. The 
volume of the solution is nearly lOOcm^ so the volume fraction occupied by Ca''', CO3 is 
(roughly) pg = 4.1 x 10^^, and our estimate of a is a ~ 7.8. The bad news is that we 
do not know /? any better than we know uj in the first place, so we cannot do primia-facie 
evaluations of A;* or the RHS of (jl.64p . 

Our policy is to use (ll.64p to obtain bounds on /? for which our model is valid. 
For instance, we have seen that the quasi-static approximation requires f3 ^ 1. In addition, 
our whole analysis is based on small super-saturation, so 77 ^ 1. This is compatible with 
(fLMj) only if 



Starting with a = 7.8, we find the LHS is unity if /? — a = 2, but if we increase the activation 
energy, (3, by Sk^T we get P — a = 5 and the LHS of p.65p is 0.05 ^ 1. Inserting this 
value of (3 into ()1.62p for k^ results in A;* ~ 1.5 x 10^. 
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Chapter 2 



The Three Eras 



2.1 Introduction 

Aggregation of identical particles (monomers) into clusters is a universal phe- 
nomenon throughout physics, chemistry and biology |KGT83j . |Kel91j . [NBOOj . |GWS+01] . 



[Isr9T], |JNpB02| , [LSgT] . |XH91j . |MG96j . |GMOJN03j . There are two classical models. The 
Becker-Doring (BD) theory [BD35j . which is based on the 'surface reactions' of adding or 
subtracting one particle (monomer) at a time from the surface of a cluster. In the Lifshitz- 
Slyozov (LS) theory of diffusion-limited growth (DLG) |LS61j a gradient of monomer concen- 
tration surrounding a cluster supports a diffusive influx of monomers into it. It is generally 
accepted that BD models the early stage of nucleation, in which 'super-critical' clusters are 
created according to the Zeldovich nucleation rate, and that DLG models the subsequent 
growth of the large clusters. In the previous chapter, we present a more general model that 
contains BD and LS as limits of small and large cluster sizes, respectively. A continuum 
limit of this expanded model yields a signaling problem for evolving the distribution of clus- 
ters in the space of their sizes. The cluster distribution satisfies an advection PDE, with the 
advection velocity specified according to DLG, and an effective boundary condition (BC) 
at zero which represents the creation of clusters at the Zeldovich rate. Both the advection 
PDE and 'creation' BC contain the excess monomer concentration, the super- saturation, 
seen far from clusters as a parameter, and the super-saturation is, in turn, a functional of 
the cluster size distribution as dictated by the overall conservation of particles. This is the 
source of strong nonlinearities in the signaling problem, nonlinearities that lead to the novel 
separation of scales behind the 'three era' mentioned in the abstract. The analysis of each 
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era consists of finding convenient variables and scales, deriving the reduced equations and 
finally solving the equations. As a result, we have the relevant physical scales of each of the 
eras. 

The chapter is structured as follows. In section 2 we formulate the signaling prob- 
lem. The governing equations and the perturbation parameter (the initial super-saturation) 
are introduced. In sections 3-5 we resolve the three distinguished limits that correspond to 
the eras: 'nucleation', 'growth' and 'coarsening'. Effective initial conditions for the growth 
era follow from matching with the 'tail' of the nucleation era, and the tail of the growth era 
similarly provides initial conditions for the coarsening era. 

Section 3 is the analysis of the nucleation era. Prom the reduced equations we 
derive a single integral equation for the super-saturation as a function of time. Prom 
the numerical solution of the integral equation we readily obtain the evolving cluster size 
distribution. In section 4, we characterize the cluster distribution during the growth era 
by ODE's for its mean and deviation. The distribution is 'narrow' because the mean size 
is large compared to the deviation, so at any given time during the growth era we have 
an ensemble of clusters all with nearly the same size. The growth slows down when the 
super-saturation is depleted to a sufficiently low level, and the mean cluster size seems to 
stabilize. However, the deviation continues to grow slowly (relative to the characteristic 
time of the growth era). This is the precursor to the coarsening era. Initially, we follow the 
slow widening analytically, by using a linearized approximation to the advection velocity 
over the (still narrow) support of the cluster size distribution. Later as this approximation 
becomes less valid due to the widening of the support, the solution to the reduced equations 
of the coarsening era is followed numerically. The numerical solution displays full-blown 
coarsening. Smaller clusters dissolve back into monomers, which are then adsorbed by the 
surviving large clusters. The numerical coarsening era solution eventually relaxes to a unique 
self-similar distribution, selected from the one parameter family of similarity solutions to 
the original equations of LS. 

2.2 The Signaling Problem 

The detailed physical-mathematical reasoning behind our reduced model of nucle- 
ation and coarsening is the subject of the previous chapter. A summary of essentials, given 
in the current section, serves as the starting point of the current chapter. 
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First, we identify the initial conditions that lead to aggregation, as opposed to 
equilibrium with negligible densities of large clusters. Equilibrium, if achievable, is charac- 
terized by densities pk of A;-clusters given by 

Pk = Pk = p\e^, k>2. (2.1) 

Here, all densities are in units of ^, where v is the molecular volume and is the binding 
energy of a /c-cluster: the energy cost to dissolve it into separated monomers. Of course, 
El = 0. For k ^ 1 we use the standard model, 

£k^kBT{ak-^aki^ . (2.2) 

The positive constants a and a are prefactors of terms proportional to the volume and 
surface area of the cluster. The total density of particles is 

oo 

P = Y.^Pk. (2.3) 

fc=i 

Substituting the equilibrium densities ()2.ip into (|2.3p gives p as as power series in monomer 
density pi. It converges for < /Oi < Ps, where ps is the critical monomer density 

Ps = e-". (2.4) 

Hence, 

oo 

fe=i 

in < pi < Ps- The total particle density, p{pi) is an increasing function. Its maximum 
value, achieved at pi = ps, is the critical particle density, pc = p{ps)- 

Since the first term of the sum in (|2.3p is pi and the rest of the terms are positive, 
it follows that pc > Ps- If we start from pure monomer with < pi < pc, we expect that 
the densities pk = Pk{t) relax to equilibrium values (j2.ip as t ^ cxd. If initially pi > pc, we 
expect aggregation and coarsening. 

During an aggregation process, the super- saturation 

_ Pi- Ps 



P: 



s 



decreases from the positive initial value ^"^^ towards zero as t oo. Positive super- 
saturation is a 'driving' parameter of aggregation. If < 0, all clusters have strong ten- 
dencies to dissolve into monomer, independently of their size. If > there is a critical 
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size 



so that sub-critical clusters with k < stiU have a strong tendency to shrink, but super- 
critical clusters with k > kc have a strong tendency to persist and grow. The essential hurdle 
to get aggregation going is a long and highly unlikely sequence of 'favorable fluctuations' 
whereby an initially sub-critical cluster eventually exceeds critical size and grows. An 
analysis along the lines of Zeldovich [Zel43j shows that the rate of creation of super-critical 
clusters per molecular volume u is asymptotic to 

jr^^e n = u;J — , (2.6) 

V DVT 

as r/ ^ O"*". Here, lo is the rate constant for a single particle to dissolve into the monomer 
bath from the surface of a cluster. We see that the creation rate is exponentially small as 
r] ^ 0^ . Nevertheless, it is not zero. Aggregation is going to happen, but it might take a 
long time. 

There are three essential ideas in our reduced model of aggregation: transport 
of clusters in the space of their size, their creation at the Zeldovich rate (|2.6p . and the 
functional dependence of super-saturation on the distribution of cluster sizes as dictated by 
overall conservation of particles. The latter is essential, because the super-saturation has a 
strong effect upon both the creation rate and the transport of clusters. 

Transport is modeled by diffusion limited growth (DLG) in which the size of any 
particular cluster is treated as a deterministic function of time, k = k{t), which satisfies the 
ODE 

k = u{k,r]) = d{r]k^ - a). (2.7) 
The prefactor d is proportional to the diffusion coefficient of monomer, 

ci=(3-167r2)^^. 

1/3 

To describe an ensemble of clusters undergoing DLG, we introduce a continuum approxi- 
mation of the cluster densities. 

Pk{t)^r{k,t) (2.8) 

for k ^ 1. Here, r{k,t) is a smooth function of its arguments. We propose that r{k,t) 
satisfies the advection PDE consistent with individual cluster sizes evolving according to 
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the ODE ((221). That is, 

dtr + dk{ur) = 0. (2.9) 

In the previous chapter, we show that DLG emerges from a full model, which contains both 
surface reaction and monomer diffusion, in the large cluster limit 

A;> A;* = (3 • levr^)!:^. (2.10) 

In all this, we are tacitly assuming a preponderance of k ^ k^.. clusters, all undergoing DLG. 
The argument for this is related to the creation of super-critical clusters, which we take up 
now. 

The creation of super-critical clusters at the Zeldovich rate |Zel43| is represented 
by an effective BC on r at A; = 0: 

1 _ 

drjkar^ne 2^, asfc^O. (2.11) 

The heuristic argument behind (j2.1ip is given in the previous chapter. Here, we make 
two notes: First, the BC ()2.1ip is shorthand for an asymptotic equality of the Zeldovich 
rate ()2.6p and the advection flux ur in an 'intermediate' range of k, much larger that 
the critical size kc, but much smaller than the characteristic size of k associated with the 
distribution r{k,t). Notice that drjki, which appears in the LHS of ()2.6p . is the dominant 
component of n in S> ^c. Second, an exponentially small creation rate as r/ — > 0"*" leads to 
an exponentially long time interval of significant nucleation. The time for a single super- 
critical cluster to grow to size A; ^ A;* is presumably much shorter. This is why we expect 
a preponderance of clusters in the DLG range of size. A; S> A;*. 

The final ingredient of the reduced aggregation model is conservation of particles, 
which implies that the super-saturation is a functional of the cluster size distribution r{k,t). 
After an initial transient, in which quasi-static densities of sub-critical clusters (A; < kc) are 
established, the conservation is expressed by 



P 



/•oo 

Pc + Psp'{Ps)v+ / kr{k,t)dk. (2.12) 
Jo 



Here, p{pi) is the equilibrium relation (j2.5p . In the RHS, pc + PsP'{pi)v approximates 
the total density of particles in sub-critical clusters in the 77 ^ 0^ limit. The integral 
approximates the density of particles in super-critical clusters. We rewrite (|2.12p as 

1 

PsP^iPs)' 



rj ^ rj^ — A / kr{k,t)dk, A 
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Here 77=,, is an effective initial value of the super-saturation, 

^ P-Pc 

PsP'ips)' 

This value is achieved after the quasi-static, sub-critical densities are established, but before 
there is a significant number of particles in super-critical clusters. 

The asymptotic analysis of the solution uses ry^, as a gauge parameter. That is, we 

define 

e = ??* 

and scale r] with e. Replacing rj with eij yields the signaling problem (SP) 

dtr + dk{vr) = 0, in k > 0, (2.13) 
V = d{khr] - a), (2.14) 

derjk^r^VLe 2^^, as A; ^ 0+, (2.15) 

poo 

7] = l-£-^A kr{k,t)dk. (2.16) 
^0 

In the following sections we determine the e — > asymptotic solution of the SP (|2.13fl2T6|) 
with IC 

r{k,0) = 0. (2.17) 

2.3 The Nucleation Era 

During the nucleation era, almost all of the supercritical clusters are created. It 
lasts until the Zeldovich creation rate decreases to a small fraction of its initial value. First, 
we find the relevant scales of the variables from dominant balances in (j2.13I[2TT6]) . The 
reduction of PDE (|2.13|) has simple characteristics that can be used to translate the SP 
into an integral equation for S7](t), the change of the super-saturation from its initial value. 
Finally, using the numerical solution for 5r](t) and the characteristics, we reconstruct the 
density of cluster sizes, r{k,t). 

2.3.1 Nucleation Scalings 

For small £ the nucleation rate (j2.15p is highly sensitive to small relative changes 
in the super-saturation. We therefore expect that a small decrease in the super-saturation 
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causes it to decrease sharply. Instead of working directly with the super-saturation, rj, we 
work with the change in super-saturation, 

5ri = r]-l. (2.18) 

We find scaling units [t], [k], [r], [6r]] of the variables t, k, r, 6r] from the dominant balances 
in equations (|2.13fl2rT6]) . Since the nucleation rate is exponentially small in e, we expect 
that a long time elapses before the super-saturation decreases enough to shut down the 
creation of additional clusters. During this long time the clusters grow to exponentially 
large sizes. Thus, the dominant term in the advection velocity in ()2.14p is k^£7] and not a. 
The dominant balance associated with (j2.14p is therefore 

[k][t]-^ = de[k]^ (2.19) 

Integrating (j2.13p with respect to k, from to oo and use of BC (j2.15p implies 

r{k,t)dk = ne . (2.20) 



d 
di 







Recalling that we expect 77 to remain near 1, we see that the dominant balance associated 
with (l2:2n]l is 

[r] [k] = Qe'^ (2.21) 



The conservation equation (j2.16p can be written in terms of 6r], 

POO 

£5r] = —A / kr{k,t)dk. 
Jo 

Therefore, the dominant balance associated with the conservation equation (j2.16p is 

e[6ri] = A [kf [r] . (2.22) 

This gives us three dominant balance equations relating the four scales of the problem. To 
find unique scalings, one more equation is needed. This last equation identifies the change 
in super-saturation 6ri that gives rise to a significant relative decrease in the Zeldovich nucle- 

3 

ation rate (j2.15p . For 6r] <^ rj, the change in the exponent — ^1^1 in (|2.15p is approximately 
-^^Sr], so the relative change in the nucleation rate, ^ is 
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Since rj is already rescaled to be near 1, in order to have an 0(1) relative change in j we 
must have 



a-" 



(2.23) 



The solutions of equations (fZlQl [2:211 {TM [2:231) for the scahng units [r] , [k] , [t] , [dr]] are 
summarized in the scaling table: 



Variable 



Scaling Table 
t 



5ri 



Unit 



,] 4 o" 



' e ^ 
i Ha- 



The scaled versions of ()2.13ff2T6|) are 

dtr + dk{ur) =0, in > 0, 

u =/c3 (1 -I — -5ri) - s, 



/ 3 / 



1 H — ^Jt] I A;3r — > exp 



v 



'2^2 



5r) = — kr dk, 
Jo 

In (|2.25p . s is the exponentially small combination of parameters 



s = 



(2.24) 
(2.25) 

(2.26) 
(2.27) 



(2.28) 



In (I2.26p . the exponent reduces to 5r] as, e ^ oo, and this 'confirms' the unit ^ of 5r] in the 
scaling table. 



2.3.2 Reduced Nucleation Kinetics 

In the limit e ^ 0, equations (j2.24I[2:271) and the initial condition (j2.17p reduce to 

(2.29) 
(2.30) 

krdk, (2.31) 



dtr + dk{k3r) =0, in A; > 0, 

kir ^e^'', as A; ^ 0+, 



6r] = — 
r{k,0) =0. 



(2.32) 
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Figure 2.1: The characteristics of the scaled 
PDE. The flux ksr is constant along the char- 
acteristics. The density of the curves matches 
the value of the solution r(k, t). 



We solve the reduced SP (|2.29fl232|) by transforming it into an integral equation for the 
change in super-saturation, Sr]{t). This integral equation is solved numerically, and the 
solution for r{k,t) is recovered from dr](t). 

The flux of super-critical clusters, j = k^r, is constant along the characteristics 



k 



. 3 

2 \ 2 

-{t -ip)j , for t > V9 



(2.33) 



of the PDE (I2.29p . The characteristics can be seen in Figure [2TT1 In the figure, the (hori- 
zontal) density of the characteristics indicates the value of the distribution of cluster sizes 
at time t. Physically, each curve in (12.33P describes the world-line of a cluster nucleated at 
time t = ip. The region below the thick line in Figure \2A\ corresponds to t < |A;2; places 
where the first nucleated clusters have yet to arrive. Hence, the value of r{k, t) there is zero. 

For a given 67]{t), the solution, r{k, t), which follows from kar being constant along 
characteristics and the BC (I2.30|) . is 



r{k,t) 



0, 



3, 2 

t>-k-., 

3, 2 

0<t < -ks. 
2 



(2.34) 



The integral equation for 5r]{t) emerges by substituting (j2.34p for r(k, t) into the conserva- 
tion identity (|23T]) . We get 



5r?(t) 



t /2 







dip. 



(2.35) 
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Notice that the variable of integration in (j2.35p has been changed from k in (j2.3ip to 
if = t- |A;2/3 as determined by (p33|) . 

2.3.3 Physical Predictions 

j(t) 




Figure 2.2: The Zeldovich flux rate j as a 
function of time. After t = 0(1), the super- 
saturation decreases slightly, and the flux gets 
turned off. 

We solve (]2.35p numerically. Figure 12.21 shows the flux as a function of time. At t 
near 1 the flux is half its original value, and that at t = 5 it has effectively vanished. A short 
discussion of the method and numerical result can be found in Appendix lA. 31 Once Sr]{t) 
is determined, we recover the distribution r{k,t) of cluster sizes from (I2.34p . Figure [2131 
displays r vs. k for a sequence of times t. 

We present collateral predictions for the total density of clusters generated dur- 
ing the nucleation era, the characteristic cluster size and the characteristic time for their 
nucleation. The (scaled) density of super-critical clusters is given by 

/•oo 

^ Pk ^ J r dk. 

Using the change of variables from k to according to (j2.33p . and the PDF (j2.29p . we 
convert the integral of r into an integral of the flux j, 

poo poo 

R= rdk= j if) dip. (2.36) 
Jo Jo 

The value of R, based on the numerical approximation to j{t), is 

R ~ 1.34. 
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Figure 2.3: The density of cluster sizes, r{k,t), 
for various values of t. 



Converting back to original physical units, the total density of clusters produced 
during the nucleation era is 



R 

— |r 



(2.37) 



where g 



2^ 



is the energy cost to create a critical cluster. Similarly, the characteristic 



cluster size [k] and characteristic time of formation [t] in original physical units are 

dvt 



[k] 
[t] 



e9 Y 



(2.38) 
(2.39) 



The nucleation asymptotics is not uniformly valid as t ^ cxd. The clusters grow by 
adsorbing monomers and the super-saturation decreases according to (j2.16p . The nucleation 
era assumption that ~ 1 loses validity. In the next section, we assume a negligible amount 
of nucleation and we follow the distribution of clusters as they evolve during the subsequent 
growth era. 



2.4 The Growth Era 



The growth era begins when the nucleation of additional clusters becomes negligi- 
ble. The continued expansion of the existing clusters depletes the super-saturation, which 
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in turn slows their further expansion. The 'tail' of the growth era is characterized by small 
super-saturation, < rj <^ 1, and a narrow, almost steady distribution of cluster sizes. 

Mathematically, the asymptotic solution to the SP during the growth era is speci- 
fied by the size of the largest clusters, and the characteristic deviation from the largest size. 
The largest cluster size and characteristic deviation satisfy simple ODE's, whose solution 
can be computed explicitly. 

2.4.1 New Variables 

We introduce new variables for the location and width of the distribution. Let 
K{t) be the size of the largest cluster as a function of time. The support of r{k,t) is 
< k < K{t). In this sense, K{t) is the 'front' of the distribution. The largest cluster, 
like any other one, grows according to iiT = u{K,r]), where u is taken from (|2.25p . It is 
convenient to introduce a 'profile' variable, x, which represents translation and scaling of 
cluster size k: 

,2.40) 

a[tj 

Heuristically, the scaling function a{t) represents the characteristic deviation of cluster size 
k from the front, K{t). Its precise meaning and determination are forthcoming. The 
distribution of clusters in x-space, denoted by q{x,t), is related to r{k,t) by 

q{x,t) = ^r{K{t) + a{t)x,t). (2.41) 

The front of the new distribution, q{x,t), is located at x = and its support is in x < 0. 
Recall that R in ()2.4ip is the integral of r{k,t) defined in (|2.36p . It is included in the 
definition of q to normalize it so that 

r-o 

qdx = l. (2.42) 



/ 



Translating SP (l2M2:27j) for r{k,t) into a SP for q{x,t) yields 

k = K^r]-s, (2.43) 

= qt + iwq)^ in X < 0, (2.44) 

-ax + (K + ax)^ri - K^-q .^x 

w = ^ — (2.45) 

a 

e^R 

r] = l ^/ {K + ax)qdx. (2.46) 

CJ-l J_K 
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Here, q and w are functions of x and t, while a, K and rj are functions of t. Notice that 
the transformed advection velocity in x-space, w, vanishes at x = regardless of a{t). This 
guarantees that the 'front' of the distribution remains at x = 0. 

We now turn to the precise determination of the characteristic deviation, or 
'width', a{t). First, we can specify its time evolution so that the x-advection velocity, 
w, has ■Wx{0,t) = 0, as well as ■w{0,t) = 0. Setting WxiO,t) = in ()2.45p gives the ODE 

a = -K-^ai]. (2.47) 
Using this choice of a, the advection velocity in (j2.45p becomes 

Assuming that a <^ K, the advection velocity w in ()2.48p is ^-^ ^^2/3 ^ fo^' ^ = 0(1)- We 
show that the asymptotic determinations of a(t), K(t) and r/(t) indeed satisfy o(t) <^ K{t) 
and ?7 <C K"^^^ during the growth era, and during the 'beginning' of the coarsening era as 
well. Hence, the x-distribution q(x, t) is nearly time- independent. 

2.4.2 Scalings (relative to nucleation era) 

The super-saturation, r], goes from being close to 1 (at the tail of the nucleation 
era) to being close to zero (at the tail of the growth era). Therefore, we do not rescale rj. 
The dimensionless signaling problem (|2.43l 12.441 12.481 12.461 12.47P is based on nucleation era 
units [k] and [t] of cluster size and time, in ()2. 381 [239]) . Hence rescaling of K and t based on 
dominant balance in the SP are relative to [k] and [t]. For instance, the dominant balance 

3 3 

of terms in (|2.46p gives as a 'relative' unit of K, so the actual unit of K is [k]-^- 
Similarly, a dominant balance between K and K'^t] in (j2.43p (neglecting the exponentially 

2 2 

small s) gives the relative unit of time, ^47^2731 s-'^d the actual unit of time is [i] ^47^273- 
The natural unit of the advection velocity, w, is length over time. In our case, 'length' is x 
which has 0(1) units, and the units of time are the units of t. Therefore, we use ^^''^J^'^''^ as 
the units of w. The ODE (I2.47P is invariant under scaling of a. But, it follows from ()2.47p 
and the leading approximation to (j2.43p . K ~ K^^^r], that a is asymptotically proportional 
to K^/^. This suggests a unit for o, -^ryf^sTs • The real justification for the unit of a comes 
from the asymptotic matching between the growth and nucleation eras, where it is shown 
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that the constant of proportionahty between scaled a and scaled K^^^ is indeed a universal, 
dimensionless constant. 

The scales are summarized in the table below. 

Scaling Table 
Variable q K rj t aw 



Unit 1 1 

Here, r is defined as 

1 2 



The scaled versions of (f2:i3l [Ml [TiSl [2:i6l [TiTD are 



. 1 



= Ksri - ST, (2.49) 
qt = -Mx, (2.50) 

V = l-K f (l + T^^)qdx, (2.52) 



„ K V K 

a = -K-lar]. (2.53) 
3 

2.4.3 Reduced Growth Kinetics 

The reduced equations that govern the dynamics during the growth era, are found 
by taking the e — > limit of ()2.49fl233]) and using the normalization (j2.42p . The first four 
equations yield 

k = K^rj, (2.54) 
qt = 0, (2.55) 
r] = l-K. (2.56) 

The advection velocity in (|2.51|) is 0(r^), and so it vanishes in the e — > limit, as in (j2.55p . 
We can therefore denote q{x,t) as Q{x) during the growth era. expected. Equation (j2.53p 
for a is invariant under the limit e — > 0. Using (|2.54p we rewrite it as 3| = This can be 
integrated, resulting in the aforementioned proportionality of a to K^^^. 

a = CK^, (2.57) 
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with the constant C to be determined by asymptotic matching with the nucleation era. We 
now determine K{t). Substituting expression (|2.56p for -q into (|2.54p gives the 

k = kI{i - K). 

The solution to (|2.4.3p . subject to K{Q) = 0, is given imphcitly bj|3 

2 

t = ^rj log (l + rjK^^ . (2.58) 
j=0 

2j + l ■ 

Here, rj are the cubic roots of —1: Vj = e 3 Figure [231 shows the solution together 
with its asymptotic behavior as t ^ and t ^ oo. 

t 




Figure 2.4: The size of the largest cluster, K{t), during 
the growth era. The two dashed lines show asymptotic 
behaviors of K{t). For small values of t, K{t) ~ (f^)^ 
and for large values of t, K{t) ~ 1 — 2.1e~*. 



2.4.4 Asymptotic Matching with the Nucleation Era 

We determine C and Q{x) by asymptotic matching of the nucleation and growth 
era solutions. This is done by examining the nucleation era solution for time that is long 
relative to the nucleation time and short relative to growth time, that is <C t ^ 1. During 
the nucleation era, the density of clusters is given by (j2.34p 

r{k,t) = k-ljl^ - hi^ , j{t) = e^^W. (2.59) 
^The trivial solution where K{t) = is ignored; we assume that K{t) ^ for t > 0. 
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Here, t and k have nucleation era scalings. Writing (j2.59p in the growth era variables using 
()2.40p and (j2.4ip . and giving t, K and a growth era scahng, we get 



a f K ax 
+ — 



tR \T 



2 Vr 



r / 



(2.60) 



The asymptotic behavior of K{t) is taken from the t ^ of the growth era, 



3 

Kit) ~ ( ^t) , for t < 1. 



In appendix IA.4I we show that substituting (j2.6ip into ()2.60p results in 

C 

q{x,t) = Q{x) = —j{-Cx) , for < t < 1. 
Jx 



(2.61) 



(2.62) 



Examining (j2.62p we see that the long-time limit of the nucleation era solution, written in 
the variables q and x, is time- independent, agreeing with (j2.55p . The asymptotic matching 
does not uniquely determine q and C: different choices of C lead to differently scaled 
distribution functions q. We are free to normalize q by choosing the constant C so that 

.0 

/ x Q{x) dx = —1. 

J —oo 

This determines a unique value for C, 

/•oo 

C = R tj{t)dt ~ 0.81. 
Jo 



The value 0.81 of C is based on the numerical solution of the integral equation (j2.35p . 
Figure [231 shows the normalized distribution profile, Q{x). 




Figure 2.5: The cluster size distribution 
q{x) = ^j{-Cx), with C ~ 0.81. 
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While the choice of constant C is arbitrary, it only affects q and not the solution 
in 'real' cluster size, k. Reverting back to r{k,t), we find that the distribution of cluster 
sizes during the growth era is given by 

r(k, t) = — (r-^^^^) , (2.63) 



which is independent of our choice of constant C. In (j2.63p . the time t is measure using 
growth era scales, t~^, and k is measured using the same scales as K, r^'^. 

2.4.5 Physical Predictions 

The characteristic time, [t] growth) of the growth era in original physical units is: 



[*W..~(;^J'|^j^j (2.64) 

2 

Recall g = ^ is the energy cost to create a critical cluster. The characteristic cluster size 
is, 

[*l«,ow.K~(j^e»)'|5. (2.65) 

Since the distribution's width a ~ r~^C is much smaller than its location, K ~ t^^, the 
clusters all have approximately the same size, [A:]growth- 

2.5 Coarsening Era 

Within the growth era approximation, the cluster distribution has a definite limit 
as t ^ oo. In this limit the distribution is concentrated in an 0{t~^) interval at A; ~ 
(recall that r = ^^—^-^), and the super-saturation vanishes. 

This is not the whole story: the growth era approximation neglects the surface 
energy component, s, of the advection velocity in (j2.49p . The remaining component of 
the advection velocity is proportional to the super-saturation, which then depletes to zero 
as the clusters grow. The full advection velocity for K, which does not neglect s, (j2.43p 
has a fixed point at a slightly smaller value of K than predicted in the growth era. A 
smaller value of K implies that the super-saturation maintains a small, yet non-zero value 
in this limit, which, in turn, implies that the width of the distribution, a, will continue to 
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grow, according to (j2.53p . with a characteristic time which sets the timescale of the the 
era. The widening of the initially narrow distribution is only a precursor to the era. Once 
the distribution is not narrow, the effective advection velocity, w in (I2.5ip is no longer 
negligible and the distribution profile changes. The distribution eventually asymptotes to 
a self-similar solution of the original LS theory. 

In this section, we scale the equations and find a simple advection PDE that 
captures the dynamics of the solution during the coarsening era. The initial conditions for 
this PDE are found by matching the t — > limit of the coarsening era solution with the 
t ^ oo limit of the growth era. Using the initial conditions found by this matching process 
we evolve the coarsening era solution using a numerical solver. 

2.5.1 New Variables and Scaling 

The coarsening era is characterized by the balance between the terms in the ad- 
vection velocity (12.250 and a small super-saturation. The distribution of cluster sizes slowly 
widens. Eventually, it fills the whole interval [0, K]. Thus, the description of the distribution 
as an advecting narrow profile is not valid here, and the subsequent choice for the 'width' 
of the distribution, a in ()2.47p is not convenient for this era. Instead of concentrating on a 
narrow distribution near the largest clusters, we consider a distribution that will be as wide 
as K. Thus in (12.43I[2l46]l we now chose a so: 

ait) = K{t). (2.66) 

Consequently, x as in (|2.40p is no longer a useful variable for the distribution. Instead we 
define 

y = x + l, q{y,t) = —r{Ky,t). (2.67) 

This is equivalent to defining y = as can be seen from (j2.40p . Using this choice for a 
and the variable y in equations (|2. 431(2^^6]) gives the following equations, written using the 
nucleation scaling. 

= qt + {wq)y, for < y < 1, (2.68) 
w = riK-l (^yl -y^ +K-^s{y-l), (2.69) 

r, = l-T'K [\qdy, (2.70) 
Jo 

K = K^ri-s. (2.71) 
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As in the growth era, the conservation equation (j2.70p imphes a scale of for 
K. The scale of the super-saturation, r] is found from the dominant balance of the terms 
in the RHS of the equation for K (12.711) . which implies ij = ts The time-scale is much 
larger than that of the growth era. Dominant balance of all three terms in (|2.7ip gives the 
characteristic time Similarly, either term in the RHS of ()2.69p implies a scale sr^ for 



ST 

w. 



The scales of the variables are summarized in the following table. 

Scaling Table 
Variable rj K t w q 



Unit ST T ST^ 1 



These scales are relative to those of the nucleation era. 

2.5.2 Reduced Coarsening Kinetics 

We rewrite equations ()2.68H277T]) using the scales above. In the e ^ limit they 



read 



= qt+{wq)y, for 0<y< 1, (2.72) 

w = K-lrj(y^-y^+K-\y-l), (2.73) 

= 1- KM, (2.74) 

k = K^7]-l. (2.75) 



Here, M is the first moment of q, 

M = [ yqdy, (2.76) 
Jo 

itself a function of t. The value of the supersaturation, rj, is implicit in this system of 
equations. In the growth era, it is given explicitly by equation ()2.56p . however, in the 
coarsening era, the scales have conspired so that r] vanishes from the corresponding equation 
(j2.74p . An explicit expression for the value of rj can be extracted from the other equations. 
The details of this derivation are in appendix IA.5[ Here, we give the results. 
During the coarsening era, the super-saturation, r] is given by 

MI Mo 

3 
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where, Mq and Mi are moments of the distribution q, 

3 

Mo= [ qdy, Mi= [ y'^qdy, (2.78) 
Jo ^ Jo 

and M is defined in (127761) . Substituting (f2:771l into the PDE ([ZTSllZTHIl for 5 and replacing 
K by jj using ([221]), we find 

= qt + {wq)y, 0<y<l, (2.79) 

w = M(^^{y-^-y)+{y-l)y (2.80) 

This PDE is nonhnear due to the dependence of the advection velocity on the moments 
of the solution. The distribution that emerges from the growth era is very narrow and a 
linear approximation of the advection velocity (j2.80p results in a PDE that can be solved 
analytically. The analytic solution widens and eventually the linear approximation loses 
validity and we turn to a numerical solution. The PDE (j2.79|. 12.801 has similarity solutions 
and the numerical solution converges to one of them as t ^ 00. 

2.5.3 Asymptotic Matching with the Growth Era 

The distribution of cluster sizes that emerges from the long time limit of the growth 
era, is concentrated in a narrow distribution near y = 1. The distribution's width is 0(t^) 
and thus, not uniformly valid as e, r ^ 0: A 'point' distribution at y = 1 is stationary under 
the action of the PDE, but a distribution with any finite width widens. To find effective 
initial conditions, which are valid as e, r — > 0, we asymptotically match the coarsening era 
solution to the long-time limit of growth era. Initially, we approximate the solution by a 
narrow distribution slowly widening due to the linearization of the PDE. Eventually, non- 
linear effects must not be ignored and this approximation loses validity. At this point we 
turn to a numerical solver to follow the continued evolution of the solution. 

The long time limit of the growth era is found from (j2.63p . Taking the limit — > 1 

we find 

r{k,t) = - 

Naively, we expect this to provide the IC for q{y,0). We pursue this venue only to change 
it slightly later. Thus, from the definition of q{x, t) in (j2.4ip we have 
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and 

K = M = Mo = Ml = r/(0) = 0. 

As we see, in y-space, this is a very narrow distribution and we approximate PDE (|2.79l - 
I2.80p by its leading order term about y = 1 

qt + {wq)y = 0, (2.82) 
w = ^{y-l). (2.83) 

The solution to this PDE subject to the IC in OrsTI) is 

This widening solution is valid as long as it it narrow, that is, as long as 

e{*-*o)/3 ^ i_ 

We see that the naive asymptotic matching leads to an initial condition that de- 
pends on r and therefore, cannot be used directly as a universal solution q{y, t) that is valid 
in the coarsening era regardless of r. However, by shifting the origin of the coarsening time, 
by to, the different functions q{y, t) for different values of r, collapse onto the same function: 

The timeshift to is dependent on r, but the resulting q is not. To summarize: we shift 
time in the coarsening era so that the asymptotic matching with the growth era happens at 
t = —to = 61ogr. This results in a universal function q{y,t) that captures the distribution 
on cluster sizes in y space during the coarsening era, and is uniformly valid for all r ^ 1. 
During an initial phase, the distribution is narrow, and an analytic approximation for q is 
given by 

At t = the distribution q so wide that the linear approximation for the PDE must not be 
trusted. We therefore, must take the solution q from some time t < when the distribution 
is still narrow, and use it as IC for a numerical solver that uses the full, non-linear PDE to 
approximate q{y, t) for larger values of t. 
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2.5.4 Numerical Solution 

For the solution of the full PDE (j2.791 12.80p we turn to a numerical solver. We 
assume that the numerical solver can accurately resolve distributions that are wider than 
a finite size 6. As we have seen, the distribution, q{y,t), starts with width (presumably 
with t'^ < 6) at t = Glogr, and widens according to (j2.84p . As the distribution q{y,t) 
widens, the linear approximation to the advection velocity loses validity. Thus, S must be 
small enough so the linear approximation to the PDE is still valid when the distribution 
has width 6. 

The solution has width 5 when the shifted time, t, satisfies e*/'^ = 5. In other 
words, we start the numerical solver with the IC in (j2.84p with 

t = Slog 5. 

The resulting solution is globally valid. Different values of r will imply a longer time 
during which the distribution is widening, prior to the numerical solution. In the resulting 
distribution q{y,t) the origin of time must be shifted back byt— >t + to = ^~ Glogr. 

r 

20 
15 
10 
5 

0.5 1 1.5 2 2.5 f 

Figure 2.6: The numerical solution at various times as found using 
clawpack. The dark curve is the initial condition, found from 
()2.84p at t = —10. The horizontal axis is k (in units of t^'^)- The 
solution continues to evolve after t = 0, but the dominant change 
is a linear stretching of the distribution. 

We acquired numerical results by using LeVeque's conservation law PDEs software 
package, clawpack |LeV02) . with the Riemann solver rpladecon (a Riemann solver for 
conservative advection). To find a numerical solution using clawpack, we used a grid 
size of 1/500 for the numerical solver, and chose a value of e~^^^^ ~ 0.0357 for 5. This 
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corresponds to starting the numerical solver at t = —10. The moments, M, Mq, and Mi 

3 

are computed by the trapezoidal rule. By comparing to numerical results obtained from a 
smaller 6 and grid spacing we estimate the error to be ~ 1%. 

To reconstruct r{k, t) we inverted the relation for t) in (j2.67p using ET, found 
from (j2.74l I2.76p . This results in a distribution in /c-space, with k and t still scaled with 
the coarsening era scales, tau~^ and — . 

O ' ST 

r{k,t + to) = RMq{kM,t), to = -6 log r. 

Figure \2M shows the solution, r(k,t), at various values oft. 

Due to our definition of the origin of coarsening time, t = 0, the interesting part of 
the solution happens at —10 < t < 2. The solution at times 61ogr < t < —10 is captured by 
the widening solution (j2.84p . At large t, the numerical solution asymptotes to a similarity 
solution of the PDE (j2.79l [XBU]) . In the next section, we review the derivation of the family 
of similarity solutions and identify among it the long-time limit of the numerical solution. 



2.6 Similarity Solutions 

The coarsening era equations (j2.76ll2T78Vl2.80p have a family of similarity solutions, 
as shown by LS. Here, we review the derivation and show that the family of similarity 
solution can be parametrized by their order of contact with zero at the largest clusters. 
This allows us to predict a priori the self-similar solution to which the numerical solution 
will converge. In the end of the section, this prediction is confirmed by comparing the 
numerical solution to the predicted similarity solution. 

To find the family of similarity solutions we look for a solution q{y,t) of the PDE 
(j2.79l I2.80p using separation of variables into 'temporal' and 'spatial' functions, c(t) and 

q{y,t)=c{t)P{y). (2.85) 
First, we find the temporal part, c(t). Recall the definition of M: 

M = [ yq{y,t)dy = c{t) [ yP{y)dy. 
Jo Jo 

Hence, for a similarity solution as in (|2.85p . M is proportional to c. We denote the constant 
of proportionality by F: 

F= yPdy. (2.86) 
Jo 
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In other words we have that, 

M = cF M = cF. (2.87) 

Another equation for M can be derived from the definition of M, the PDE (j2.79|) and 
integration by parts, 

M= [ yqt = - [ y{wq)ydy= [ wqdy. (2.88) 
JO Jo Jo 

From the two equations (|2.87p and (|2.88p for M, we derive an ODE for c: 

cF = M^{l-n). 



Here, /i, is 



Thus, c satisfies the ODE 



Mo loPdy 



Ml JlyPdy 



c = -Fc'^ifj.-l). (2.89) 
Notice that < /x < 1 and that fi is time- independent, hence the solution of ODE (j2.89p is 

for some tg- 

Now that we have c{t), we can find the spatial part of the similarity solution, 
P{y). Substituting the ([2:85]) into the advection PDE p7l9\i . and using the ODE ([2:89]) for 
c, yields an ODE for P: 

P(c^F{l-f,)+cWy)+PyCW = 5i ^ _Mi:ii^l±^. (2.91) 

The explicit solution of ([2.9ip . subject to P(0) = 1, is 

/x(2-iy-|)-2 



P{x) = exp / ^ dx . (2.92) 

/U(2/3 -y) + (y-l) J 

The product of c{t) and P{y) gives a two-dimensional family of solutions parametrized by 
fi and tg. The parametrization by tg in ([2.90p reflects time-independence of the PDE 1)2.79] 
I2.80p . Different values of select different similarity solutions for fi G [0, 1]. However, we 
only allow values | < ^ < |. Other values of /x generate P that explodes at one of the 
end-points, and we dismiss such solutions as non-physical. The parametrization by tg in 
([2.90p reflects time-independence of the PDE ([2.791 12.80p . We use this symmetry during the 
asymptotic matching to the numerical solution. Figure \27l\ shows a few similarity solutions 
for different values of /i. 
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Figure 2.7: The plot of the similarity solution for 
various values of fi. The solutions are normal- 
ized so that Mq = 1. The dark line has ^ = |. 



2.6.1 Order of Contact at y = 1 



Here, we find which of the similarity solution profiles in (j2.92p is the limit of the 
numerical solution in the coarsening era. For this we study the order of at y = 1 contact 
of the various similarity solutions. The similarity solutions are originally parametrized by 
fi. It it more instructive to parameterize them by their order of contact with zero. At 
the tail-end of the coarsening era, the numerical solution is discontinuous at y = 1. This 
corresponds to an order of contact 0. Thus, we find the similarity solution with order of 
contact zero, and confirm that the numerical solution converges to it. 

The order of contact is the power p so that 

P(y)~6(l-yf asy^l", 
for some constant 6 7^ 0. A little algebra shows that 

Pyiy 



p = lim (InP)y(l — y) = lim 
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P 



(2.93) 



Substituting (j2.9ip into (|2.93p and using I'Hopital's rule we find that the order of contact 
of the similarity solution P{y) is 

f = ^- (2.94) 



2// -3 

Thus, we expect the numerical solution of the coarsening era to converge to the similarity 
solution with p = 0, hence 



6 
5' 



Equation (j2.9ip is integrable for this value of /x, resulting in 

125 exp (-2^1 (coth-^ (^21) - tanh"! (^^f^ 



P{y) 



X 



2/3 



(2.95) 



(2.96) 
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Here, P{y) is normalized so that 



Pdy = l, 



yPdy ~ 0.632573, 



'0 Jo 
The dark, discontinuous hne in Figure [2771 shows P with fi 



1 1 5 
.yiPdy = -. 
6 



5- 



(2.97) 



2.6.2 Asymptotic Matching with Coarsening Era 



1 . 5 



0.5 




Figure 2.8: The numerical solution q(y,t), 
scaled so that M = 1, at different times. Start- 
ing on the right at i = —10 in a narrow distri- 
bution, and converging to the similarity solution 
with /U = I (Dark line) at t = 4. 

As can be seen in Figure [2T8| the numerical solution of the coarsening era converges 
to the similarity solution with /i = | as t ^ oo. The temporal behavior of the similarity 



solution given by (j2.90p contains an undetermined constant, tg- By asymptotic matching of 
the numerical solution to the similarity solution we now determine this last undetermined 
constant. 

From (p:871) for M and (pIOO]) for c{t) it follows that the first moment of the 
similarity solution decays like 

1 



M(t) 



(2.98) 



(/i-l)(t-t,)' 

for some ts- Hence, a numerical solution that converges to the similarity solution with 
^ = |, must satisfy 

j^ = t-ts, ast^oo. (2.99) 
From the numerical solution we calculate M{t) and plot t^tv. By fitting a line of 



slope 1 to the linear part of the plot of jj we find that tg 



M{t) ■ 

8.5. The results are plotted 



in figure E 
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Figure 2.9: The match of to t — tg. The 



approximate numerical value of tg is -8.5. 



Figure 12.91 shows that the long-term behavior of the numerical solution's first 
moment matches the predicted decay as t — > oo. The same can be said about the profile of 
the numerical solution. Figure 12.81 show the convergence of the (scaled) numerical solution 



numerical solution should be close to the similarity solution. Indeed, at t = 5 the numerical 



2.6.3 Physical Predictions 

Our asymptotic solution to the LS PDF predicts the emergence of a distribution 
of clusters which is discontinuous at the large clusters. The predicted similarity solution, 
has a specified timeshift tg. This does not include the timeshift Iq which is needed for the 
initially narrow distribution to widen prior to the numerical solution. The scaling gives us 
specific physical scales for the size of the clusters, [k]c, and the time of their formation, [t]c- 



to the similarity solution with /i = |. The figure seems to indicate that at t > 5 the 



solution differs from the A* = | similarity solution by less than 1%. 




(2.100) 



(2.101) 



2.7 Conclusions 



The three eras, nucleation, growth, and coarsening emerge from the asymptotic 
analysis of the aggregation model which combines elements from BD and LS models. At 
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the small super-saturation limit the three eras are separated by increasingly longer times 
scales. The solution starts with IC r{k,0) = 0, corresponding to pure monomer. Each 
era provides effective initial conditions for the next one by means of asymptotic matching, 
and the coarsening era selects the discontinuous similarity solution as the global longOtime 
limit of the aggregation process. The selection of the similarity solution is accompanied by 
physical scales and a specific shift ts in (j2.90p . both were previously unknown. Due to the 
vast difference in scales between the different eras, we do not expect a single experiment to 
resolve the complete nucleation process, from nucleation to coarsening. However, different 
experiment could be set-up to examine each of the different eras by adjusting the physical 
parameters. 

Future work is possible in several aspects of this solution: 

1. Post-coarsening effect of the "ignition" and diffusion 

The ignition phase, during which the nucleation rises from zero to the rate specified 
by the Zeldovich formula, takes a small, yet finite time. This implies that the front 
is not sharp, but rather transitions smoothly to zero. In addition, in deriving the 
advection PDE a diffusion term was ignored because it is asymptotically small at all 
3 eras. Both effects imply that the order of contact with zero at the front is not 0, 
and that there is be another scale of k in which the distribution transitions smoothly 
to zero near k ^ K. Therefore, at a time-scale even larger than that of the coarsening 
era, this 'contact layer' can widen and possibly change the selection of the long-term 
similarity solution. Whether or not this is physically relevant depends on the timescale 
and in which this happens. 

2. A successor to the Zeldovich nucleation rate 

The Zeldovich nucleation rate is calculated using a mean-field for the monomer dis- 
tribution around a cluster. This approach does not take into account the role of 
fiuctuation in the monomer density. Aggregation is a discrete process, and nucleation 
a rare-event; the nucleation rate might be found by taking this into account. 

3. Corroboration with experiments 

The scales predicted in the present chapter have not been checked against data from 
experiment. The main difficulty in doing this is that in the experiment the large 
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clusters tend to fall out of the solution due to gravity. There have been some exper- 
iments in micro-gravity and it would be informative to compare our results to the 
experiments. 

4. Extension to two dimensional models 

The present thesis models space and clusters as 3-dimensional. This assumption sets 
the exponents governing the advection velocity of the PDE. Two-dimensional crystals 
can be grown on silicon wafers and other flat surfaces. Other experiments show growth 
of three dimensional crystals on silicon wafers. Crystals growing on a two-dimensional 
substrate will form large clusters without gravity interfering and will allow for long 
experiments. Reformulating the current model for other dimensions and comparing 
with experimental data could be fruitful. 

5. Inhomogeneous far-field supersaturation 

Our model assumes that the far-field supersaturation is constant. This has been shown 
to be an unstable solution. A model which allows for the far-field supersaturation to 
vary in space will present a more realistic solution. 
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Appendix A 

Appendix 



A.l Deterministic Cluster Growth According to BD 

Here we determine the domain of k in which the average change in the cluster size 
{5k) is a good approximation to the actual change in cluster size. BD assumes that the 
cluster grows and shrinks by means of discrete, independent Poisson processes. The cluster 
gains monomers at a rate CkPi via the adsorption process, and via the emission process it 
loses monomers with rate d^. Therefore, the average change in size over a short interval 5t 
is 

{5k) = {k{t + 5t) - k{t)) = {ckPi - dk)5t, 
provided that 5t is small enough so that 

\ckPi — dk\ 5t <^ k. 

Since the two controlling process are assumed to be independent Poisson processes, 
the mean square deviation of k from the average is 

{{5k-{5k)f) = {ckPi + dk)5t 

The process can be considered deterministic if the deviation is much smaller than the 
expected change, so 

V^cfcpi +dkV5i < I (cfcpi - dk)5t\ . (A.2) 
Combining inequalities (jA.ip and ()A.2p we get a condition on 5t 

CkPi + dk<. \{ckPi - dk)f 5t <C \ckPi - dk\ k. 
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We can find 6t which satisfies this condition if 



CkPi +dk <. \ckPi - dk\k. 



(A.3) 



This does not hold for every value of k. In particular, we notice that for k = kc the RHS 
is zero. Hence, the k{t) that can be approximated as deterministic must be bounded away 
from kc- 

For A; ^ 1 in ()A.3P we insert the approximations (jl.5l II. 6|) for dk and , and then 
examine its asymptotic form as 7/ — > 0+. The result is 



2 + r,|l-( A 



< r] 



Jxif 



k. 



(A.4) 



For clusters with k ^ k^, this condition is satisfied trivially for t] <^ 1. For k on 
the order of kc, we need to dig a little deeper. 

Since we are looking at the 77 — > limit, we look at the leading order terms of 

foil 

2 

1 



rjk 



kc\ 3 



For k near kc we expand the RHS around k = kc- Keeping the leading order terms of the 
Taylor series of the RHS, we find an explicit condition for the evolution of clusters to be 
treated deterministically: k must be bounded away from kc according to 

- < A; - /cc . 

As r/ — > 0, the excluded domain grows. On the other hand, its size relative to kc ^ 



shrinks to zero: 



a" 



Therefore, we see that for small super-saturation the evolution of clusters can be treated 
deterministically, for all but a vanishingly small domain around kc- 



A. 2 Effective Super-Saturation, 77, After Ignition 

It has previously been shown that an initial condition of pure monomer leads to a 
transient 'ignition' phase in which quasi-static densities of sub-critical clusters are created 
|NBC05j . Since these sub-critical clusters are made of monomers, after the ignition transient 
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the monomer density will be lower than the original value, prior to the ignition. We compute 
the value rj^: of super-saturation after the formation of the k < kc quasi-static densities, but 
before significant depletion by super-critical clusters. 

The starting point if the sum (jl.39p truncated to k < kc because the contribution 
from super-critical clusters is insignificant: 

p = ^kpk. 

l<k<kc 

Since the sub-critical cluster densities are quasi-static, we approximate them by the equi- 
librium densities (|1.46p . and the corresponding approximation of the sum is 

p = ^/cpje"^. 

l<k<kc 

Recalling the definition of the super-saturation, (jl.Op . pi = /9s(l + ^*), and taking the 
two-term expansion as 77* — > 0'*', we find 

l<fc<fcc 

The sum on the RHS can be obtained by differentiating the series representation ()1.47p of 
p{pi), setting pi = ps, multiplying by ps, and then truncating to k < kc. The derivative 
series converges at pi = ps, so for rj ^ (and hence for kc 00), the sum is asymptotic to 
Psp'{ps). Hence, 

P ~ Pc + 'n*PsP{Ps), 

and so, 

1 P- Pc 
P'ips) ps 

A. 3 Numerical Solution of Integral Equation 

The nucleation era solution required an integral equation to be solved. The scaled 

equation 

determines the change in the supersaturation with time. We found the approximate solution 
of this equation on a discrete set of point separated by a fixed interval h. The value at each 
point was calculated from the trapezoidal rule approximation to the integral of the solution 



dip 



(A.5) 
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up to time t + h. This is not implicit, since the t — (p term in the integrand vanishes at 
If = t. The order of accuracy of the method is found to be | by numerical experiment. This 
fractional order is most likely due to the cusp in the integrand which adds an error of /l3/2 
to the integral. 

The value R=1.343 was found by using 1000 points in the interval (0, 4). The error 
in R is estimated to be 0.001, while the error in the flux itself is estimated to be 10~^. The 
errors were obtained by comparing to the results of the same numerical calculation with 4 
times as many points. 



A. 4 Asymptotic Matching Between Nucleation and Growth 
Eras 

Here, we examine q{x, t) in the overlap domain <^ t <^ 1 between nucleation 
and growth given by (j2.60p 

a ( K ax\ .1 t 3 f K axy^ 



For small t, K{t) has asymptotic behavior given by (j2.6ip . and a{t) is a simple function of 
K, as given by 

3 

a = CK\. (A.7) 
Selectively substituting the expressions in (IXtI) for K{t) and a{t) into (IXBl) yields, 
, , CkI K-\ ( r'^axY^ ft kI( Zt'^qxY' 



Cancellations and additional use of (lA.Tjl give 

Cf T'^axY^^ ft f / 3t^CxY' ' 

Therefore, as r — > 0, t fixed 

q{x,t)=Q{x) ~ J^ji-Cx). 
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A. 5 Finding 77, the Details 

Equations (j2.72H2776]) describe the evolution of the particle size distribution im- 
plicitly. The the supersaturation, rj, is needed for the advection velocity (|2.74p . but it has 
no explicitly defining equation. 

To find rj, we derive another equation for K from (j2.74p and the PDE (|2.72l I2.73P 

1 ■ M 

The value of M can be found from the definition of M ()2.76p , the advection equation (|2.72p , 
and integration by parts 

^ = y yQtdy = J wqdy =^ K = j wqdy. (A.9) 

Comparing the two expressions for K in ()2.75p and ()A.9P gives and expression for rj, 

, = Mi(l-£^), ,A.10) 

which is implicit because w ()2.73p depends on rj. To untangle the dependence on rj, we 
write Jq wq dy using rj and various moments of q. 



2 / / 1 



wqdy = K'i7]l [y'i - y] qdy + K ^ j {y-l)qdy 
JO ^ ^ JO 

2 



K"3r?(Mi -M) +K"i(M-Mo). (A.ll) 



Where M, Mq and Mi are moments of the distribution q, 

3 

M = I yqdy, Mq = / qdy, Mi = [ y^iqdy. 
Jo Jo ^ Jo 

Substituting U^TMh and (IXTTl) in dXlO]) for ry results in 



Mar? Ml - M + M(M - Mq) 

77 = Mt 11 — ^ ^— 

' I M2 



which can be massaged into a much simpler expression for 77: 

3 



Hence, we end up with the following advection PDE to solve: 

= qt + {wq)y, 
w = M 



Mo 
M 



^ (y3 -y) + (y- 1) j 



